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Abstract 

In this paper, we consider the Dirichlet problem of inhomogeneous incompress- 
ible nematic liquid crystal equations in bounded smooth domains of two or three 
dimensions. We prove the global existence and uniqueness of strong solutions with 
initial data being of small norm but allowed to have vacuum. More precisely, for 
two dimensional case, we only require that the basic energy Hy^uoll^ + [|Vdo 11^2 
is small, while for three dimensional case, we ask for the smallness of the production 
of the basic energy and the quantity ||Vuo||^ 2 + ||V 2 <io||^2- Our efforts mainly cen- 
ter on the establishment of the time independent a priori estimate on local strong 
solutions. Taking advantage of such a priori estimate, we extend the local strong 
solution to the whole time, obtaining the global strong solution. 

Keywords: existence and uniqueness; global strong solutions; liquid crystal. 



1 Introduction 

We consider the following hydrodynamic system modeling the flow of nematic liquid 
crystal materials 

p t + div{pu) = 0, (1.1) 
p(u t + (u ■ V)u) - uAu + Vp= -Adiv(Vd Vd), (1.2) 
divw = 0, (1.3) 
d t + (u- V)d = 1 (Ad+\Vd\ 2 d) (1.4) 

in Q x (0,oo), where Q is a bounded domain with smooth boundary in Mr (N = 2,3). 
Here u : X (0, oo) —> M N represents the velocity field of the flow, d : ft x (0, oo) — > S 2 , the 
unit sphere in IR 3 , represents the macroscopic molecular orientation of the liquid crystal 
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material, p : Q x (0, oo) — > [0, oo) and p : Q x (0, oo) — )■ R are scalar functions, respectively, 
denoting the density of the fluid and the pressure arising from the usual assumption 
of incompressibility divu = 0. The positive constants u, A and 7 represent viscosity, 
the competition between kinetic energy and potential energy, and microscopic elastic 
relaxation time or the Dehorah number for the molecular orientation field, respectively. 
The symbol Vd © V<i, which exhibits the property of the anisotropy of the material, 
denotes the N x N matrix whose (i,j)-th entry is given by did ■ djd for 1 < i, j < N. 
Noticing that 



System ( ll.ll) - (11.4l) is a simplified version of the Ericksen-Leslie model, which reduces 
to the Ossen-Frank model in the static case, for the hydrodynamics of nematic liquid 
crystals developed by Ericksen [TJ, [2] and Leslie [3] in the 1960's. Both the full Ericksen- 
Leslie model and the simplified version are the macroscopic continuum description of the 
time evolution of the materials under the influence of both the flow velocity field u and 
the microscopic orientation configurations d of rod-like liquid crystals. A brief account 
of the Ericksen-Leslie theory and the derivations of several approximate systems can be 
found in the appendix of [I]. For more details of physics, we refer the readers to the two 
books of Gennes-Prost [5] and Chandrasekhar [6] . Though the above system is a simplified 
version of the full Ericksen-Leslie system, it still remains the most important mathematical 
structures as well as most of the essential difficulties of the original Ericksen-Leslie system. 

In the homogeneous case, i.e. p = C, Lin-Lin [4j [7] initiated the mathematical analysis 
of f ll.2p - fll.4p in the 1990's. More precisely, they considered in [I] the Leslie system of 
variable length, that is replacing |V<i| 2 <i by the Ginzburg-Landau type approximation 
term d to relax the nonlinear constraint \d\ = 1, and proved the existence of global 
weak solutions in dimension two or three. They also obtain the unique existence of 
global classical solutions in dimension two or in dimension three with v large enough. 
Furthermore, they proved in [7j the partial regularity theorem for suitable weak solutions, 
similar to the classical theorem by Caffarelli-Kohn-Nirenberg jS] for the Navier-Stokes 
equation. A preliminary analysis of the asymptotic behavior of global classical solutions 
was also given in [1] . More precise asymptotic behavior of classical solutions can be found 
in Wu [9], in particular, he provided an estimate on the convergence rate in dimension 
two. However, as pointed out in [U [7], both the estimates and arguments in these two 
papers depend on e, and it's a challenging problem to study the convergence as e tends 
to zero. The two dimensional case is comparatively easier, in fact Hong [10J obtains the 
convergence as e goes to zero up to the first singular time. Such convergence problem in 
dimension three is still open. Alternatively, one can establish the existence of global weak 
solutions directly to the system fll.2p -f jl~4|) but for the Ginzburg-Landau approximate 




one can rewrite equation 11.21 as 




XAd ■ Vd. 
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system. Recently, Lin-Lin- Wang [TT] proved the global existence of weak solution to the 
system (jl.2p — (11.4p in dimension two, and obtained the regularity and asymptotic behavior 
of the weak solutions they established. The uniqueness of such weak solution was later 
proven in [12]. For three dimensional case, the local or global existence of weak solutions 
is still unclear in the present. 

In the non-homogeneous case, i.e. the density dependent case, the global existence 
of weak solutions to the system (ll.ip ~ fll.4p with |V<i| 2 <i being replaced by 1 2 d, the 
Ginzburg-Laudan type approximation term, is established in [T2J [T4] and [T3] for each 
e > 0. They cannot get the uniform estimates with respect to s > 0, and therefore cannot 
take the limitation e — » 0. It's also a challenging problem to study the convergence as e 
tends to zero for the non-homogeneous case. If the initial data gains more regularities, 
one can expect to prove the existence of more regular solutions. In fact, Wen and Ding 
[j~6] obtain the local existence and uniqueness of the strong solutions to system fll.ip - fll.4p 
with initial density being allowed to have vacuum. If the initial data is small or satisfies 
some geometric condition, one can obtain the global existence results: global existence of 
strong solutions in three dimensions with small initial data are obtained by Li and Wang 
in [T7] for constant density case, Li and Wang in [18] for nonconstant but positive density 
case, and Ding, Huang and Xia in [19] for nonnegative density case; global existence of 
strong and weak solutions in two dimensions with large initial data is obtained by Li [2U] 
under the condition that the third component of the initial direction filed is away from 
zero. 

In the present paper, we consider the global existence of strong solutions to the liquid 
crystal equations. More precisely, we establish the global existence of strong solutions to 
the non-homogeneous system fll.ip - fll.4l) . coupled with the following initial and boundary 
conditions: 

(p,u,d)\ t=0 = (p ,u ,d ), with |d | = 1, divu = 0, and u \ dn = 0, (1.6) 




where d$ is a given unit constant vector and po(x) a given nonnegative function being 
allowed to vanish on some subset of Q. Compared with the approximation problem, the 
term |V<i| 2 <i in fll.4l) brings us some new difficulties, for example, one can not obtain the 
a priori L 2 estimates on Ad from the basic energy identity. System fll.ip - fll.4p can be 
viewed as Navier-Stokes equations coupling the heat flow of harmonic maps. Since the 
strong solution of a harmonic map can be blow up in finite time [21] , one cannot expect 
that (jl.ip - (jl.7p have a global strong solution with general initial data. Therefore, we 
consider the case that the initial data is of small norm but the initial density po is allowed 
to have vacuum. 

Throughout this paper, for any 1 < p < 00, we denote by ||w|| p the LP norm of u for 
any u G L P (Q). Using this notation, we can state the main result of this paper as follows: 



Theorem 1.1. Let Q be a bounded smooth domain in IR (N = 2,3). Assume that 
po G H\n) n L°°(n), < p (x) < p in SI, u G H 2 (tt) n ^(O) with dWu -Oiufl, 
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do G H 3 (Q) with \do\ = 1 and do = d^ on dfl, for some positive constant p and a constant 
unit vector dy, and the following compatible condition is valid 

-vAuq - Vp - Adiv(Vd Vd ) = y/p 9o 

in n for (p ,g ) e H\tt) x L 2 (VL). 

Then there is a positive constant e depending only on v, \, 7, p and Q, such that if 

\\\fpo~u \\l + ||Vd ||a < £ o for N = 2, 

or (\\Vfru \\l + \\Vd \\ 2 2 )(\\Vu \\l + ||V 2 rfo||2) < ^0 for N = 3, 

then initial boundary value problem ( ti. 1\) - (UT\) has a unique global strong solution (p, u,p, d) 
satisfying 

p G L™ ([0, 00); H\n)) n L°°(0, 00; L°°(fi)), G L~([0, 00); L 2 (fi)), 

u g l°°(o, 00; # 2 (ft) n i/o( fi )) n l2 (°> °°; ^ 2,6 (^)), 

^ G L 2 (0, 00; H^tt)), ^pu t G L°°(0, 00; L 2 (n)), 
p G L°°(0, 00; H\n)) n L 2 (0, T; W 1 ' 6 ^)), 

d G W 2 4 ' 2 (g) nL°°(0,oo;if 3 (fi)), rf t G L°°(0, 00; ^(fi)), |d| = 1, 
where Q = Q x (0, 00). 

We now comment on the analysis of this paper. Since the local existence of strong 
solutions to system f ll.ll) - fll.7l) has been proven in [16], to establish the global existence 
result, we only need to extend the local solution to the global one. For this aim, recalling 
that vacuum is allowed in our paper, we need to establish some a priori estimates on 
local strong solutions, which is independent of the existence time interval and the lower 
bound of the density. The first key estimate of this paper is the estimate on Ei(t) (see 
Lemma 13.21 for the definition), which controls the L°°(0,T; H 1 ) norm of the velocity u 
and the L°°(0,T; H 2 ) norm of the direction field d. Via energy estimates, we obtain a 
polynomial inequality of Ei(t) with small coefficients (see Lemma [3721 for the detail), with 
an additional term involving J Q * ||u||oo|| VwH^s if N = 2. This additional term results from 
the assumption imposed on the initial data that only the initial basic energy is small for 
N = 2, and it disappears if we impose the same assumption for iV = 2 to that for N = 3. 
On account of such polynomial inequality on Ei(t), one can use continuity argument 
to derive the bound of Ei(t) if iV = 3, while for N = 2, we can employ a logarithmic 
type Sobolev inequality to obtain the bound of Exit). As long as we obtain the estimate 
on Ei(t), the next step is to do higher order estimates, i.e. the estimates on E 2 (t), 
which controls the L°°(0, T; H 2 ) norm of the velocity u and the L°°(0, T; if 3 ) norm of the 
direction fields d. Similar to the situation encountered before, the arguments are different 
for N = 2 and those for N = 3, and the cause is still the different assumptions imposed on 
the initial data for N = 2 and N = 3. In fact, if we use the same approach to that used for 
iV = 3 to the deal with the case iV = 2, we will encounter a term J* \\Vd\l\\u t \\lds, which 
can not be controlled in terms of Ei(t) under the assumption that only the basic energy 
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of the initial data is small. After obtaining the higher order estimates on local strong 
solutions, we use the standard approximation approach to establish the global existence 
of strong solutions, and thus finish the proof. 

The rest of this paper is arranged as follows: in Section [2} we state some preliminary 
lemmas which will be used in the next two sections; in Section [3j we do the a priori 
estimates on the local strong solutions, including both the basic energy estimates and 
the estimates on higher derivatives of u and d independent of the existence time interval 
and the lower bound of density; in Section HI taking advantage of the a priori estimates 
established in Section [31 we prove the global existence and uniqueness of strong solutions 
by using the standard approximation approach. 

Since the exact values of v, A and 7 don't play a role, we henceforth assume 

v = X = 7 = 1 

thoughout this paper. We denote 

Co = WVPoUqWI + HVdoll! 
the basic energy of the initial data in the rest of this paper. 

2 Preliminaries 

In this section, we give some useful lemmas which will be used in the rest of this paper. 
The following result is quite standard (as a matter of fact, it's a straightforward 
generalization of the one presented in |22j). 

Lemma 2.1. Let Q be a Lipschitz domain ofM> N and v £ L 1 (0,T; Lip) be a solenoidal 
vector-field such that v ■ n — on dQ, where n denotes the outward normal vector on dfl. 
Let p £ W 1,q {Vt) with q G [1, 00]. Then equation 

{p t + div(pv) =0 in O, 
p\t=o = Po inVL 

has a unique solution in L°°(0, T; W 1,0 °(n)) D C([0, T]; ni< r<00 W 1 ''"(Q)) if q = 00 and in 
C([0,T]; W^iQ)) ifl<q<oo. 

Besides, the following estimate holds true 

\\p(t)\\w^(n) < eti^ v{T) ^ dT \\p \\w^n) 
for any t £ [0, T]. If in addition p belongs to L p for some p £ [1, 00] , then 

\\p{t)\\ P = IIPo|| P 

for all t £ [0, T]. Finally, if p (x) > 5 in Q for some positive constant S, then p(x,t) > 5 
for all (x,t) £ x [0,T]. 
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We will frequently use the following two lemmas, which state the elliptic L q estimates 
on Laplace and Stokes equations. 

Lemma 2.2. (See JE3S) Let Vt be a C k+2 bounded domain in R N , N > 2, k > is an 
integer. Then we have 

\\u\\ W k+2, q ^ < C(\\ Au\\ W k,q(fy + ||w||i9(n) + ||<7||w fc + a '«(«)) 

for any u G W k+2,q {Vt) and g G W k+2,q {Vt) with u\qq = g, 1 < q < oo, where C is a 
positive constant depending only on N,q and Q. 

Lemma 2.3. (See J '24]) Let tt be a bounded domain in R N , N > 2, of class C m+2 , 
m > 0. For any 

f G W m ' q (n), if G w m+2 - 1/q ' q (dn), 

1 < q < oo, with compatible condition j dn if ■ ndS = 0, there exists one and only one pair 
u,p such that 

(i) vE W m+2 > q (tt) and p G W m+1 > q (n)/R, 

(ii) v,p verify the Stokes equation 

—Au + Vp = / in Q, 
div u = in Q, 
u\dn = 

In addition, this solution obeys the inequality 

\\ v \\w m + 2 'i{n) + ||p|| w m + 1 -i{n)/R < C(||/||w m >«(n) + |Mlw m + 2 - 1 /</.9(afi))) 

where C is a positive constant depending only on N,m,q and Q. 
We also need the following local existence result. 

Lemma 2.4. (See JIE$) Under the conditions stated in Theorem \l.l\ there is a constant 
T* > 0, such that for any T < T* system ( ti. 1\) - (UT\) has unique solution (p,u,p,d) 
satisfying 

p e L°°(0, T; H\Q)) H L°°(0, T; L°° {£!)), p t G L°°(0, T; L 2 (Q)), 

u G L°°(0,T;# 2 (ft) nH^(n)) nL 2 (o,T;iy 2 - 6 (fi)), 
u t G L 2 (0,T;H^)), ^pu t G L°°(0, T; L 2 (Q)), 
p G L°°(0, T; H\Q)) H L 2 (0, T; W 1,6 (fi)), 

d G W^(Q T ) n L°°(0,T; H 3 (Q)), d t G L°°(0,T; H^Q)), \d\ = 1, 
where Q T = Q x (0,T). 



3 A priori estimates 

In this section, we concern on the energy estimates on strong solutions. Let T > 
and (p,u,p,d) be a strong solution to (ll.ll) - (ll.7p in Q x [0,T) stated in Lemma [23 By 
Lemma [2. 1[ we have 

< p(x,t) < p 

on Qt- Recalling that d\on = d^, then Sobolev embedding theorem implies 

\\V m d\\l = \\V m (d-d*)\\l<\\d-d* \\ 2 Hm+1 < \\V m+1 d\\ 2 2 

for any integer m > 1 if iV = 3. We will frequently used these facts without any further 
mentions later. 

Lemma 3.1. Let (p,u,p,d) be a strong solution to U. l\) - (Ll\) in Q T . Then we have 
the following energy estimates 

sup (||VHl2+HV« + 2 [ T (\\Vu\\ 2 2 +\\Ad\\ 2 2 )dt<Co + 2 [ T \\Vd\\idt, 

0<t<T Jo Jo 

and 

sup (Hv^ulll + \\Vd\\l) + 2 / \\Vu\\ 2 2 dt < C . 

0<t<T Jo 

Proof. Multiplying (11. 2p by u, using (11. ip and integration by parts, we obtain 



d f p\u 
di 



( — — ^—dx + [ | Vu\ 2 dx — — f (Ad ■ Vd) ■ udx 
Jn 2 J n J n 

u ■ V)d ■ Addx = - Ad - (Ad + \Vd\ 2 d - d t )dx, 



from which we get 



~y- / {p\u\ 2 + \Vd\ 2 )dx+ [ (\Vu\ 2 + \Ad\ 2 )dx 
2 dt J n J n 

= - I \Vd\ 2 Ad- ddx= [ \Vd\*dx, 
Jn Jn 

here we have used the fact that Ad ■ d = —\Vd\ 2 , which is guaranteed by = A|<i| 2 
2Ad - d + 2\Vd\ 2 . Hence 



T 

2 i II A j||2\ 



sup (\\^pu\\ 2 2 + \\Vd\\l) + 2 / (\\Vu\\ 2 2 + \\Ad\\ 2 )dt 

0<t<T Jo 

<2 f \\Vd\\\dt + (||v£ uo||2 + HVdollI) < Co + 2 f \\Vd\\\dt. 
Jo Jo 
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Since |Vd| 2 = —Ad ■ d and \d\ = 1, it follows || V<i||| < ||Arf|||, and thus we deduce from 
the above inequality that 



sup (||VHl2+l|V« + 2 / (\\Vu\\ 2 2 +\\Ad\\ 2 2 )dt<C + 2 / ||Ad||^t, 

0<t<T Jo Jo 

which implies 

sup (UVHI2 + l|Vd||l) + 2 / \\Vu\\ 2 2 dt < C Q . 

0<t<T Jo 

The proof is complete. □ 

Lemma 3.2. Let (p,u,p,d) be a strong solution to U.1\) - [T77\ ) in Qt, and set 

r-t 

7 2 r||2\ , / /II /-. ||2 , ||V72„,||2 1 ||V7„||2 , ||V73j||2\ 



E 1 (t)= BUp(||V«||^ + ||V^)+ / (||^ll2 + l|V^|| 2 ^ + ||V^ + ||V^)^- 

0<s<t Jo 

Then we have the following 
(%)IfN = 2, there holds 

Exit) < C(||Vu ||l + \\V 2 d \\l) + CCoE^t) + CC 1/2 [ ||u||^||V«||^ + C [ \\V 2 d\\jds. 

Jo Jo 

(ii) If N = 3, there holds 

Exit) < C(||Vt*o||l + \\V 2 d \\l) + CCoExit) 2 + C f || V 2 d\\ 2 2 ds. 

Jo 

Proof. Multiplying ( 11.21) by Ut, integration by parts and using Young inequality, we 

get 

d f \^ u \ dx+ f p \ Ut \i dx= I (yd®VdS7u t - p{u-V)u- u t )dx 
dt Jq 2 J n Jq 

=4- f VdQVdVudx - [ ((Vdt Vd + Vd Vdt)Vu + p(u-V)u-u t )dx, 
dt Jn Jn 



and thus 



It I (^~2~ ~ Vrf ° WVm ) dx + J P\ u t\ 2<ix 

<2 / (\Vd\\Vd t \\Vu\+ p\u-Vu\\u t \)dx 
Jn 

<- / p\u t \ 2 dx + C / (\Vd\\Vd t \\Vu\ + p\u-Vu\ 2 )dx, 
2 Jn Jn 

which gives 

f (|Vw| 2 -2Vd0 VdVu)dx+ [ p\u t \ 2 dx<C [ (|Vd||Vd t ||Vu| + p\u ■ Vu\ 2 )dx. 
dt Jn Jn Jn 

(3.1) 



Applying elliptic estimates of Stokes equations to (11. 2p yields 

\\^ 2 u\\l + \NP\\1 < C{\\ P {u t + u ■ Vu)\\l + || Vd\V 2 d\\\l). (3.2) 
Combining (13.11) with (13.21) . there holds 

sup ||Vw||l+ / (IIVMII2 + \\V 2 u\\l+ \\Vp\\l)ds 

0<s<t Jo 

<C|| V^olll + C sup |||Vd| 2 Vu||i + e f \\Vd t \\lds 

0<s<t Jo 



+ / (\Wd\ 2 \Wu\ 2 + p\u\ 2 \Wu\ 2 + \Wd\ 2 \W 2 d\ 2 )dxds. (3.3) 
Jo Jn 

Taking the operator A to both sides of equation (jl.4p and then multiplying the resulting 
equation by Ad, we deduce 



d f \Ad 

dt 



2 

2. 



dx — / A dAddx 
Jn 

[ (| Vd| 2 | Ad| 2 + 2Vd : VAddAd + 2| V 2 d\ 2 dAd 
Jn 

- (Au ■ V)dAd - 2{VuidiV)dAd)dx 



<e / (|VAd| 2 + \Au\ 2 )dx + C / (|W| 2 |V 2 d| 2 + \Wu\\W 2 d\ A )dx. (3.4) 
Jn Jn 

In the above we have used the fact that d ■ Ad = — \Vd\ 2 guaranteed by \d\ = 1. Note 
that Ad\gn = |V<i| 2 c/|^ guaranteed by equation (jl.4p and the boundary condition (ll.7p . 
Integration by parts gives 

- f A 2 dAddx = - f Ad-?- AddS + f \VAd\ 2 dx 
Jn Jan on J n 

= [ \\7d\ 2 d-?-AddS + [ \VAd\ 2 dx 
Jon dn J n 

= I (\Vd\ 2 ^-(dAd) - \Vd\ 2 Ad-^d\ dS+ I \VAd\ 2 dx 

= - I [\Vd\ 2 ^\Vd\ 2 + \Vd\ 2 Ad^d\dS+ I \VAd\ 2 dx, 

which, combined with (13 .4p . it follows from the trace inequality that 

f \*f-dx+ [ \VAd\ 2 dx 
dt Jn 2 J n 

<e [ (|VAd| 2 + \Au\ 2 )dx + C [ \Vd\ 3 \V 2 d\dS 
Jn Jan 
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+ C I (|W| 2 |V 2 d| 2 + \Vu\\V 2 d\ 2 )dx 
Jn 

<e [ (\VAd\ 2 + \Au\ 2 )dx + C\\\Vd\ 3 \V 2 d\\\ w ^ {Q) 
Jn 

+ C I {\Vd\ 2 \V 2 d\ 2 + \Vu\\V 2 d\ 2 )dx 
Jn 

<e [ (\VAd\ 2 + \Au\ 2 )dx + C [ (\Vd\ 3 \V 2 d\ + \Vd\ 2 \V 2 d\ 2 + \Vd\ 3 \V 3 d\)dx 
Jn Jn 

+ C f (\Vd\ 2 \V 2 d\ 2 + \Vu\\V 2 d\ 2 )dx 
Jn 

<e I (|VAd| 2 + \Au\ 2 )dx + C I (|W| 2 |V 2 d| 2 + |W| 3 |V 3 d| 



+ \Vu\\V 2 d\ 2 + \Vd\ 3 \V 2 d\)dx 

<e [ (|V 3 rf| 2 + \Au\ 2 )dx + C [ {\Vd\ 2 \V 2 d\ 2 + \Vd\< 
Jn Jn 

+ \Vu\\V 2 d\ 2 + \V 2 d\ 2 )dx, 

and thus 



sup ||V 2 d|| 2 + f \\V 3 d\\ 2 2 dt 

0<t<T Jo 

<C\\V 2 d \\ 2 2 + e [ T \\V 2 u\\ 2 2 dt + C [ T (\Vd\ 2 \V 2 d\ 
Jo Jo 



+ |Vu||V 2 d| 2 + |W| 6 + \V 2 d\ 2 )dxdt. 
Combining (13.31) with the above inequality, and using Lemma I3.1[ we obtain 

sup (|| win + ||v 2 rf|| 2 ) + A||vWl2 + 11^112 + H v pII' + W^ d \\l) ds 

0<s<t Jq 

<C(\\Vuo\\l + \\V 2 do\\ 2 2 ) + C sup |||VtZ| 2 V«||i + C f [ (|W| 2 |V 2 d| 2 + |W| 2 |Vu| 2 

o<s<t Jo Jn 

+ pM 2 |Vm| 2 + |Vu||V 2 d| 2 + |W| 6 + \V 2 d\ 2 )dxds 

=C{\\Vuo\\l + \\V 2 do\\l) + ^Ph + C f \\V 2 d\\ 2 2 ds. (3.5) 

i=i Jo 

We estimates the terms on the right hand side of the above inequality as follows: If iV = 3, 
then it follows from Soblev embedding inequality and Lemma 13.11 that 



h <C sup ||V«|| 2 ||Vd||J /2 ||Vd||J /2 < CC 1/4 sup ||Vu|| 2 ||V 2 d||5 /2 

0<s<* 0<s<t 
7 2 7 1 1 2 i „, iiv7„.ii4 



<e sup \\V 2 d\\ 2 2 + CC sup \\Vu\\% 

0<s<t 0<s<t 



h <C I ||Vci|| 2 ||Vci|| 6 ||V 2 rf||^ S < CCl' 2 sup ||V 2 d|| 2 / \\V 3 d\\ 2 2 ds, 

JO 0<s<t Jo 

h<C f \\Vd\\ 2 \\VdU\Vu\\lds<CCl 12 sup ||V 2 d|| 2 f \\V 2 u\\jds, 

Jo 0<s<t Jo 

h<C [ \\y/pu\\ 2 \\u\\ 6 \\Vu\\lds < CC 1/2 sup ||Vu|| 2 / \\V 2 u\\ 2 2 ds, 

Jo 0<s<t Jo 

h<C f \\Vu\\ 2 \\V 2 d\\lds<C f \\Vu\\ 2 \\V 2 d\\y 2 \\VH\\Tds 
Jo Jo 

<e [ \\V 3 d\\ 2 2 ds + C [ \\Vu\\ 4 2 \\V 2 d\\ 2 2 ds 
Jo Jo 



< niV'dWlds + CCo sup \\Vu\\l\\V 2 d\\l 

0<s<t 
t-T 

72 j||6j„ / n I IIV72j||4||v7Jl|2||v73 j\\2. 



h <C / \\W 2 d\\lds < C / ||V^||*||V^||V^||^ 
Jo Jo 

<CC sup ||V 2 d|| 2 / \\V 3 d\\ 2 2 ds, 

0<s<t Jo 

if N = 2, then 

h <C sup \\Vu\\ 2 \\Vd\\\<C sup ||V«|| 2 ||W|| 2 ||V 2 d|| 2 

0<s<t 0<s<t 

<CCl 12 sup ||V M || 2 ||V 2 rf|| 2 , 

0<s<t 

h<c f \\vd\\l\\v 2 d\\lds<c f \\vd\\ 2 \\v 2 d\\ 2 2 \\v 3 d\\ 2 ds 

Jo Jo 



ft 



h<c 



[ \\Vd\\ 2 2 \\V 3 d\\ 2 2 ds<CCo [ \\V 3 d\\ 2 2 ds, 
Jo Jo 

\\Vd\\l\\Vu\\ 2 4 ds < C f ||Vrf|| 2 ||V 2 rf|| 2 ||V M || 2 ||V 2 M || 2 rf S 
Jo 

<e f \\V 2 u\\ 2 2 ds + C f \\Vd\\l\\V 2 d\\l\\Vu\\ 2 2 ds 
Jo Jo 

<e [ \\V 2 u\\ 2 2 ds + CCl sup ||V 2 rf|| 2 , 

Jo 0<s<t 

h<C f \\^pu\\ 2 \\u\U\Vu\\lds < CC'J 2 f \\u\U\Vu\\ 2 \\\V 2 u\\ 2 ds 
Jo Jo 

<e [ \\V 2 u\\jds + CC I ||u||^||V 2 u||^s, 
Jo Jo 

f \\Vu\\ 2 \\V 2 d\\lds <C f \\Vu\\ 2 \\V 2 d\\ 2 \\V 3 d\\ 2 ds 
Jo Jo 



h<C 
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<e [ \\V 3 d\\ 2 2 ds + C f \\Vu\\l\\V 2 d\\lds 
Jo Jo 

<e [ ||V 3 rf||^s + CC sup ||V 2 rf||2, 

JO 0<s<t 

h<C f \\Vd\\ 2 2 \\V 2 df 2 ds < C [ \\Vd\\i\\V 3 d\\ 2 2 ds <CC 2 [ \\V 3 d\\ 2 2 ds. 
Jo Jo Jo 

Substituting the above inequalities into (13.51) and setting 

E 1 (t) = sup (|| Vug + \\V 2 d\\l) + [\\\^pu t g + \\V 2 u\\ 2 2 + \\Vp\\l + \\V 3 d\\ 2 2 )ds, 

0<s<t Jo 

we obtain that, if N = 3 then 

Ei(t) <C(||V«o||2 + II V 2 rf ||l) + CCoE^t) 2 + CC^E^tf' 2 + C [ \\V 2 d\\ 2 2 ds 

Jo 

<C(||V«o||2 + II V 2 rf ||l) + CC E 1 (t) 2 + ^(t) +C f \\V 2 d\\ 2 2 ds, 

2 Jo 

and thus 

Exit) < C(||Vuo||! + II V 2 d„||i) + CC Q E l {t) 2 + C f \\V 2 d\\ 2 2 ds, 

Jo 

and if N = 2, then 

Ei(t) < C(||Vw ||l + \\V 2 d \\ 2 2 ) + CCoE^t) + CC^ 2 [ \\u\\ 2 J\Vu\\ 2 2 dt + C [ \\V 2 d\\ 2 2 ds. 

Jo Jo 

The proof is complete. □ 

Before continuing the energy estimates, we cite the following Sobolev inequality of 
logarithmic type, which will be used in Lemma 13.41 

Lemma 3.3. (See 125^ ) Assume Q is a bounded smooth domain in M 2 and f G 
L 2 (s,t;H 1 (n)) n L 2 (0,T; with some q > 2 and < s < t < oo. Then it 

holds that 

||/|U 2 (s,t;L°°(n)) < C(l + ||/||L2( s ,t;//i(^))( ln+ II / 1 U a (s,t; V^La (f2)) ) )> 

with some constant C depending only on q and Q, and independent of s,t. 
Now, we state and prove the following lemma. 

Lemma 3.4. Let (p,u,p,d) be a strong solution to U.l\) - fiT77]) in Qx- Let Ei(t) be 
the function defined in Lemma \3.S[ Then there is a positive constant e depending only 
on p and Q, such that 
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(i) IfN = 2, then 

E 1 (t)<C[l + C (\\Vu \\l + \\V 2 d \\l) 2 }, 

provided \\y/p u \\l + ||Vd ||a < e o- 

(ii) If N = 3, then 

E 1 (t)<2C(\\Vu \\l+\\V 2 d \\l), 
provided (\\^p u \\l + ||Vd |||)(|| Vuolli + \\V 2 d \\ 2 2 ) < e . 

Proof. We first consider (i). It follows from Lemma l37fl and Ladyzhenskaya inequality 
that 

/ \\V 2 d\\ 2 2 ds<CC + C [ \\Vd\\ids < CC + C sup \\Vd\\ 2 2 [ ||V 2 d|| 2 ds 

JO JO 0<s<t Jo 

<CC + CC f \\V 2 d\\ 2 2 ds, 
Jo 

and thus J* 1 1 \7 2 g? 1 1 < CC , provided e is small. On account of this fact, using Poincare 
inequality, it follows from Lemma 13.21 that 

Bk{t) < C(\\Vu \\i + || V 2 d \\ 2 2 ) + CCl' 2 f \\u\U\Vu\\ldt. 

Jo 

By Lemma 13.31 it follows from Poincare inequality and Lemma 13.11 that 

J Wuf^dsKC ( 1 + J o \\Vu\\ 2 2 dt) H 1 + J o (l|V«||l+ \\V 2 u\\ 2 2 )ds] 

<CIn(l + f \\V 2 u\\ 2 2 ds). 
Jo 

On account of this inequality, it follows from Gronwall inequality that 

E x {t) <C{\\Vuo\\l + l|V 2 rfo||^)C 1/2 f \\u\\ldse cc « 12 $>\\^s 

Jo 



<C(||V M o||^ + ||V 2 rfo||^)C 1/2 ln(l+ [ ||V 2 n|| 2 rf S )(l+ f \\ V 2 u\\ 2 2 ds) cc o 

Jo Jo 

<CC 1/2 (||Vu || 2 + ||V 2 d ||^)(l + f \\^ 2 u\\lds) 1/2 



<\Ex{t) + CC^(\\Vu \\ 2 2 + || V 2 c?o|||) + CC (\\Vu g + \\V 2 d 0ll2) , 
provided Eq is small enough, and thus 

Ex(t) < C[l + C7o(||V«o||5 + HV^oll 2 .) 2 ], t e (0,T). 
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Now, we prove (ii). Using Gagliado-Nirenberg, it follows from Lemma [3. II that 

f \\V 2 d\\ 2 2 ds <CC + C [ \\Vd\\lds<CC + C [ \\ W|| 2 || V 2 d\\ 3 2 ds 
o Jo Jo 

<CC + C [ ||W|| 2 ||V 2 d|| 2 ||W|| 2 ||V 3 d|| 2 d S 
Jo 

<CC + l [ \\V 2 d\\ 2 2 ds + C f \\Vd\\i\\V 3 d\\ 2 2 ds, 
z Jo Jo 



and thus 



\\V 2 d\\ 2 2 ds <CC + CC sup ||Vd||^ / ||V 3 d||lds 

0<s<t Jo 



<Ci\\Vuo\\ 2 +\\V 2 d \\ 2 ) + CC sup \\V 2 d\\i / WV'dMds 

0<s<t Jo 

<C(\\Vuo\\ 2 2 + \\V 2 d \\ 2 2 ) + CCoE^t) 2 , 

in the above we have used the fact that ||w|| 2 + || Vd||| < C(|| Vm|| 2 + || V 2 ^) 2 ,) guaranteed by 
Poincare inequality and the boundary condition ( II. 7p . On account of the above inequality, 
by Lemma |3.2| there is a positive constant > 1 depending only on p and Q, such that 

Exit) < a(\\Vu \\ 2 2 + || V 2 rf || 2 ) + aCoE^t) 2 , 

which implies 



P m < 1- >/l - 4CgC (|| Vriolll + l|V 2 rfo||l) 
Ei{t) < 

or 



E (f) > 1 + VI - 4C7gCo(||V^ ||l + ||V 2 ^o||l) 

Note that Ei(t) is a nondecreasing and continuous function on [0,T]. Set Bo — sc^' t^ en 
it has 

4C 2 Co(||Vn || 2 + ||V 2 d ||^)<l/2. 

One can easily check that 

*(o) = iiwoiii + iiv*ib < ^m™ >. 

iO + t/0 

Consequently, the continuity of Exit) implies that 



Exit) ^ 1 ~ VI ~ 4 CgCo(||Vtz ||| + ||V 2 do||l) 



2C C 

l-(l-4^C (|lWo°| ]l+||V 2 do|]l)) 
IC^Co 

the proof is complete. □ 



< vx ^T;°" 2 " °" 2 " = 2a(||v« ||| + iiv^oii 2 .), 
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Lemma 3.5. Let (p,u,p,d) be a strong solution to U.1\) - (T77\ ) in Qt- Set 
E 2 (t) = sup (\\^pu t \\i + ||V 2 M || 2 + || Vp||| + || V 3 d|| 2 ) 

0<s<t 

+ [\\\Vu t \\l + \\d u \\l + \\V%\\l)ds. 
Jo 

Then there is a positive constant e depending only on p and Q, such that the following 
hold 

(i) IfN = 2, then 

E 2 (t) < C(l + ||s || 2 + \\V 2 u \\ 2 2 + ||V 2 cZ Q || 2 H1 ) 8 e (1+ ll Wo| ^ +l|v2do ^ )2 

for all t G (0, T), provided Ha/PoMoU! + ll^ 2( ^o||2 < £ o- 

(ii) IfN = 3, then 

E2(t)<C[\\g \\l + (l + \\Vu \\l + \\V d \\ 2 H2 r} 

for all t e (0,T), provided (H^/wMl + II V 2 rf |ll)(|| V?i || 2 + \\V 2 d \\ 2 ) < e . 

Proof. Differentiating the momentum equation (jl.2p with respect to t, multiplying 
the resulting equation by u t and integrating over fl yields 

= / (Vd t QVd + VdQ Vd t ) : Vu t dx - / p(u t ■ V)u ■ u t dx 
Jn Jn 

+ / div(pu)(u t + (u • V)w) • u t dx. 
Jn 

It follows 

div(pu)(u t + (u ■ V)w) • u t dx - 

Jn 

<C / {p\u\\Vu t \\u t \+ p\u\\Vu\ 2 \u t \+ p\u\ 2 \V 2 u\\u t \+ p\u\ 2 \Vu\\Vu t \)dx. 
Jn 

Combining the above two inequalities yields 

^ f p\u t \ 2 dx+ [ \Vu t \ 2 dx 
at Jn Jn 

<C [ (\Vd\ 2 \Vd t \ 2 + p\u\ 2 \u t \ 2 + p\u\\Vu\ 2 \u t \ 
Jn 

+ pM 2 |V 2 w|M +pM 4 |Vw| 2 + pM 2 |Vw|)cfe. (3.6) 

Sine the estimates on d are different for the case N = 2 and the case N = 3, we prove 
(i) and (ii) separately. 



J jpn ■ V(M 2 ) + pu ■ V[(u ■ V«) ■ u t ]}dx 
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(i) The case that N = 2. Taking A on both sides of equation (jl.4p . there holds 

Ad t - A 2 d = \Vd\ 2 Ad + 2{Vd : V Ad)d + 2\V 2 d\ 2 d 

- (u ■ V) Ad - 2(V«i • diV)d - (Aw • V)d. 

Squaring both sides of the above equation and integrating over Q yields 

-2 I Ad t A 2 ddx+ [ (\Ad t \ 2 + \A 2 d\ 2 )dx 
Jn Jn 

<C [ (\u\ 2 \VAd\ 2 + \Vu\ 2 \V 2 d\ 2 + \Vd\ 2 \Au\ 2 



Jn 

+ \Vd\ 4 \Ad\ 2 + |W| 2 |VAd| 2 + \V 2 d\ A )dx. (3.7) 

By equation (11.41) and the boundary value condition (11.71) . there holds 

&d\an = -\Vd\ 2 d\ dn = -dl\Vd\ 2 \ dn . 

By the aid of this boundary value condition, it follow from elliptic estimates and La- 
dyzhenskaya inequality that 

J \V 4 d\ 2 dx = J \\V\d-d* )\\ 2 2 dx < C\\A(d- d*)\\ 2 H2 

=C\\Ad\\ 2 m < C(\\A 2 d\\ 2 2 + |||Vd| 2 ||| 2 ) < C(\\A 2 d\\ 2 2 + |||W| 2 || 2 + ||V 2 |W| 2 || 2 ) 
<C(||A 2 rf|| 2 +||Vrf||l+|||Vrf|V 3 rf|| 2 + |||V 2 rf| 2 || 2 ) 
<C(\\ A 2 d\\ 2 2 + \\Vd\\i + || W|| 2 || V 3 d|| 2 + || V 2 d\\t) 

<C(\\A 2 d\\ 2 2 + ||Vrf|| 2 ||V 2 rf|| 2 + ||Vd|| 2 ||V 2 d|| 2 ||V 3 d|| 2 ||V 4 d|| 2 + ||V 2 d|| 2 ||V 3 d|| 2 ) 

<i||V 4 d|| 2 + C(||A 2 rf|| 2 +||Vd||^||V 2 d||^ 
+ ||Vrf|| 2 ||V 2 rf|| 2 ||V 3 rf|| 2 + ||V 2 rf|| 2 ||V 3 C /|| 2 ), 

and thus it follows from Lemma [3.41 that 

|| V 4 rf|| 2 < C|| A 2 rf|| 2 + C(l + HWolll + || V 2 d || 2 ) 2 (l + || V 3 rf|| 2 ). (3.8) 

Note that the following boundary condition holds true 

Ad t \ 9 n = -(\Vd\ 2 d) t \ dn . 

Integration by parts, it follows 

-2 / Ad t A 2 ddx = - 2 [ -^-Ad-Ad t dS + 2 [ VAd t VAddx 



n Jan un Jn 

d 

~dt 



4- [ \VAd\ 2 dx + 2 [ -^-Ad(\Vd\ 2 d) t dS. (3.9) 
dt J n J 9n dn 
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Recalling the trace inequality 

< C\\f\\ w i,^ n) , 

there holds 

f ^-Ad(\Vd\ 2 d) t dS = [ -^-Ad(\Vd\ 2 d* ) t dS 
Jan on J dn dn 

=2 f ^-AdWd : VdtdtdS = 2 f -^-Ad ■ dWd : Wd t dS 
Jan dn J m dn 

= f (^-(Ad ■ d) - Ad^-d) 2Vd : VdtdS 
J an \9n dn J 

= - f \^-\Vd\ 2 + Ad^-d\2Vd:Vd t dS 
Jan \dn dn J 

> - C I \V 2 d\\Vd\ 2 \Vd t \dS > -C\\\V 2 d\\Vd\ 2 \Vd t \\\ m ,i {n) 

Jan 

> - C I {\V 2 d\\Vd\ 2 \Vd t \ + \V*d\\Vd\ 2 \Vd t \ 

Jn 

+ \v 2 d\ 2 \vd\\vd t \ + \v 2 d\\vd\ 2 \v 2 d t \)dx 

> - e\\V 2 d t \\ 2 2 - C [ (|W| 2 |V 2 d||W t | + |W| 2 |V 3 d||W 4 | 

Jn 

+ \Vd\\V 2 d\ 2 \Vd t \ + \Vd\ 4 \V 2 d\ 2 )dx. 
Combining the above inequality with (I3.7p - fl3.9p and using Lemma [3 A\ we obtain 
sup ||V 3 d||^ + (||V 2 rf t || 2 +||V 4 rf|| 2 ) 

0<s<t 

<C[||VArf || 2 + 1 + (||V 2 u || 2 + \\V 2 d \\ 2 2 ) 4 ] + C [\\u\ 2 \V 3 d\ 2 

Jo 

+ |W| 2 (|V 2 u| 2 + |V 3 rf| 2 + \V 2 d\\Vd t \ + |V 3 d||W 4 |) 

+ |Vw| 2 |V 2 cf + |V 2 rf| 4 + |W| 4 |V 2 d| 2 + \Vd\\V 2 d\\Vd t \]dxds. (3.10) 

By Lemma [3.11 it follows from Ladyzhenskaya inequality that 

f \\V 2 d\\ 2 2 ds < C [ \\Ad\\ 2 2 ds < CC Q + C [ \\Vd\\ 4 ds 
Jo Jo Jo 

<CC + C sup ||W|| 2 / \\V 2 d\\ 2 2 ds < CCo + CC [ \\V 2 d\\ 2 2 ds, 

0<s<t Jo Jo 



and thus 



\\V 2 d\\ 2 ds < CC , 
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provided Eq is small enough. On account of this inequality, using Ladyzhenskaya and 
Gagliado-Nirenberg inequality, by Lemma I3.1[ we estimate the terms on the right hand 
side of inequality (I3.10p as follows 

I 1= f [ \u\ 2 \V 3 d\ 2 dxds < [ \\u\\l\\V 3 d\\lds 
Jo Jn Jo 

<C [ \\Vu\\ 2 2 \\V 3 d\\ 2 \\V 4 d\\ 2 ds < e [ WV'dWlds + CE^tf, 
Jo Jo 

I 2 = I j |Vrf| 2 (|V 3 ci| 2 + \V 2 d\\Vd t \ + \V 3 d\\Vd t \)dxds 
Jo Jn 

<c r||vrf|ll(liv 3 d||l + ||vd t ||l + ||v 2 0^ 

Jo 

<C f ||Vd|| 2 ||V 2 d|| 2 (||V 3 rf|| 2 ||V 4 rf|| 2 + IIVdtHallV^Ha + || V 2 d|| 2 || V 3 d\\ 2 )ds 
Jo 

< £ [\\\V A d\\l + \\V 2 d t \\l)ds + C f || V 2 rf|| 2 (|| V 3 d|| 2 + || Vd.ll 2 + || V 2 d\\ 2 2 )ds 



<e A||v 4 rf|| 2 + ||v 2 ^|| 2 )^ + c( J E; 1 (t) + J E; 1 (t) 2 ) + J E; 1 (t) t \\vd t \\ids, 

Jo Jo 

\Vd\ 2 \V 2 u\ 2 dxds<C [ ||Vd||^||V 2 «||^s 




o Jn 



<C f \\Vd\\ 2 \\\\V 3 d\\ 2 \\V 2 u\\ 2 2 ds <e sup \\V 3 d\\ 2 2 + C ( t \\V 2 u\\ 2 2 ds 

Jo 0<s<t \Jo 



<e sup \\V 6 d\\i, + CE x {tY, 

0<s<t 



= f I (\Vu\ 2 \V 2 d\ 2 + \V 2 d\ i )dxds < f (\\Vu\\l\\V 2 d\\l+ \\V 2 d\\l)dx 
Jo Jn Jo 

<C r(||V M || 2 ||V 2 rf|| 2 ||V 2 M || 2 ||V 3 rf|| 2 + ||V 2 rf|| 2 ||V 3 rf|| 2 )^ < CE^t) 2 , 
Jo 

= 11 |W| 4 |V 2 cf cfaxfe < / \\Vd\\t\\V 2 d\\lds 
Jo Jn Jo 

< f ||Vd|| 2 ||V 2 d||2||V 2 d|| 2 ||V 3 d||l<fo < CE 1 (t) 3 , 
Jo 

= [ f \Vd\\V 2 d\\Vd t \dxds< f ||Vrf||4||V 2 rf|| 4 ||Vrfi|| 2 ^ 
Jo Jn Jo 



Jo Jn 

, 1/2 ,, _ 2 1/2 1| ^2 Ti, 1/2 ,, _ 3 1/2 



<c / iivdii^liv^iiriiv^ll^liv^n^iiv^ii^ 

Jo 

<C f {\\V 2 d\\ 4 2 + \\V 3 d\\l + \\Vd t \\l)d S <CE 1 {t) + C f \\Vd t \\ 2 2 ds. 
Jo Jo 
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Substituting all the above inequalities into (I3.10P yields 
sup ||V 3 ^ + (||V 2 ^ + \\V 4 d\\l) 

0<s<t 

<C[\\VAd \\ 2 2 + 1 + (|| V 2 u \\ 2 2 + || V 2 doHl) 4 ] + CE^t)) 3 

+ C(l+E x (t)) [ \\Vd t \\lds. (3.11) 
J o 

By equation (ll.4p . using Ladyzhenskaya inequality and Sobolev inequality, recalling that 

/ \\V 2 d\\jds < CC , 
Jo 

it follows from Lemma 13.41 and Lemma 13.11 that 

\\Vd t \\lds <C [\\\V 3 d\\ 2 2 +\\V{u-Vd-\Vd\ 2 d)\\l)ds 
Jo 

<C [\\\V 3 d\\l + \\u\\l\\V 2 d\\l + || V«||2|| W[|| + ||Vd||jj + \\Vd\\l\\V 2 d\\l)ds 



<C J (||V 3 d|| 2 + ||V M || 2 ||V 2 d|| 2 ||V 3 d|| 2 + ||V M || 2 ||V 2 M || 2 ||W|| 2 ||V 2 d|| 2 

+ ||Vrf|| 2 ||V 2 rf||^+||Vrf|| 2 ||V 2 rf|| 2 ||V 3 rf|| 2 )rf S 
<c£(||V 3 d||^ + ||V 2 n|| 2 + ||V M ||^||V 2 rf|| 2 

+ II Vn|| 2 || V 2 d\\\ + || V 2 rf|| 2 + || V 2 d\\t)ds < C(l + E^t) 2 ), (3.12) 
from which, recalling (13.111) and using Lemma 13.41 we arrive 

sup ||V 3 d||^+ Al|V 2 ^|| 2 + ||V 4 d|| 2 )d S 

0<s<t Jo 

<C[||VArf ||l + (1 + II V 2 u \\ 2 2 + ||V 2 d ||^) 6 ]. (3.13) 

By flEU), ^ has 



sup ||^ t || 2 + / ||V^|| 



0<s<t 



<\\VP~ou t m\l+ I [ (\Vd\ 2 \Vd t \ 2 + p\u\ 2 \u t \ 2 + p\u\ 2 \V 2 u\\u t \ 
Jo Jn 

+ p\Vu\\u t \ 2 + p\u\\Vu\ 2 \u t \ + p\u\ 4 \Vu\ 2 )dxds. (3.14) 

By Ladyzhenskaya inequality, Sobolev embedding inequality and recalling (13. 12ft . we can 
estimate the terms on the right hand side of the above inequality as follows 



h< [ [ {\Vd\ 2 \Vd t \ 2 + p\u\ 2 \u t \ 2 + p\u\ 2 \V 2 u\\u t \)dxds 
Jo Jn 



<C I (\\Vd\\l\\Vd t \\l + \\u\\l\\^pu t \\ 2 \\u t \U + IM| 8 ||V 2 u|| 2 |M| 4 )^ 
Jo 

<C [ (||Vd|| 2 ||V 2 d|| 2 ||Wi|| 2 ||V 2 ^ 2 + \\Vu\\l\\^pu t \\ 2 \\Vu t h 



+ ||Vw|| 2 ||V 2 w|| 2 ||V^|| 2 )ds 

<^\l|V 2 rfdl2+l|V«dl2)^ + ^£(l|V42l|V 2 rf||2l|V^|l2 

+ \\Vu\\i\\^u t \\l + \\Vu\\l\\VM 2 2)ds 
<e [\\\V 2 ck\\l + \\Vu t \\l)ds + C(l + E 1 (tf), 




|2. 



p|Vm||m 4 | dxds < C / HVwimiVpwtlhllwtlUrfs 
'o Jn Jo 

<C / ||Vu||J /2 ||V 2 u|| 2 /2 llV^tlbllVutllads 
Jo 

<e / ||Vn t || 2 d S + C / ||V«|| 2 ||V 2 «|| 2 ||VMIl2^ 
Jo Jo 

/"V^ll2^ + C A||V W || 2 +||V 2 «|| 2 2 )||VP«t|l2^ 

Jo Jo 

<£ / \\Vu t \\ 2 2 ds + C [ IIV^II^II^II^ + C^W 2 , 
Jo Jo 

h — / p\u\ | Vu| 2 |wt|(is < C I | ] 1 1 4 1 1 Vu|| 2 ||-U f || 4G?S 

Jo Jn Jo 

<C / ||VM|| 2 ||V 2 M|| 2 ||Vt^|| 2 rfs < £ / ||VM 4 || 2 rfs + C / ||VM|| 2 ||V 2 M|| 2 ds 

Jo Jo Jo 

<e [ WVutgds + CE^t) 3 , 
Jo 

h = [ [ p\u\ A \Vu\ 2 dxds < C [ ||«||^||Vw|| 2 ds 
Jo Jn Jo 

<C [ \\Vu\\t\\Vu\\ 2 \\V 2 u\\ 2 ds < C [ ||Vw|| 2 1 (||V«|| 2 + ||V 2 w|| 2 )cfe 
Jo Jo 

<C{(1 + E 1 {tf). 

Substituting these inequalities into (13.141) . using the inequality (13 . 13f) and recalling 
compatible condition, it follows from Lemma [3.41 that 



sup ||V^*ll2+ / || Vw t || jds 



0<s<t 



<C[\\g \\l + HVAdbll! + (1 + \\V 2 u \\ 2 2 + || V 2 rf ||l) 6 ] + C f ||V 2 u|| 2 ||V^ t |l2^, 

Jo 
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from which, by Gronwall inequality, we deduce 
sup / ||VtAt|||rfs 

0<s<t Jo 

<C[\\VAd \\l + (1 + || V 2 u || 2 + || V 2 d \\ 2 2 ) 6 ]e c ^ l|v2u|l ^ s f \\V 2 u\\\ds 

<C(1 + \\9o\\l + W^oWl + l|V 2 rf ||^) 6 (l + HV^olll + llV^ll^e^Hli+llv^lll) 2 
<C(1 + \\9o\\l + \W 2 u \\l + ||V 2 rfo||^) 8 e (1+ " v ^"^" v2d0 "^ 2 . (3.15) 
Elliptic estimates to Stokes equations give 

liv 2 «|| 2 + ||v P || 2 <c(\\^u t \\i + |Mv M ||| 2 + ||vrf|v 2 rf||| 2 ) 

^(llv/p^ll^+ll^ll^llVnll^+llVrfll^llV^II 2 ) 

<C(\\V~m\\t + || V M || 2 ||V M || 2 || V 2 M || 2 + II W|| 2 ||V 2 rf|| 2 ||V 3 d|| 2 ) 

<e(|| V 2 «|| 2 + || V 3 rf|| 2 ) + C(||^|| 2 + + ||V 2 rf||^), 

which, combined with (13 . 131) . together with (13. 15ft . we obtain 

sup (H^ll 2 + ||V 2 M || 2 + ||Vp|| 2 + || V 3 rf|| 2 ) + A|| V^|| 2 + || V 2 ^|| 2 + IIV^H 2 )^ 

0<s<t Jo 

<C(1 + \\go\\l + II V 2 ^o||l + ||V 2 rf ||^) 8 e (1+l|Wol|2+llv2dol|2)2 . 

(ii) The case that N = 3. Differentiate equation (jl.4p with respect to t, then it has 

d tt - Ad t = \Vd\ 2 d t - (u ■ V)d t + 2(Vd : Vd t )d - (u t ■ V)d. 

Square both sides of this equation and integration by parts, using Young inequality, we 
have 

j I \Vd t \ 2 dx+ I (\d tt \ 2 + \Ad t \ 2 )dx 

<C [ (|W| 4 K| 2 + M 2 |W*| 2 + |W| 2 |W t | 2 + |W| 2 K| 2 )cfe. (3.16) 
Jn 

Combining this inequality with (13. 6p . we obtain 

sup (\\VpMI + ||Vrf t || 2 ) + A||Vn t || 2 + \\d tt \\ 2 2 + \\Ad t \\ 2 2 )ds 

0<s<t Jo 

<(IIV^«*(0)ll2 + l|Vd t (0)|| 2 ) + C f [ (plnpKp + plnHV^H^I + lnHVrf^+lVrfpKp 

Jo Jn 

+ I Vrf| 2 | Vrftl 2 + pM 4 |Vw| 2 + \Vd\ 4 \d t \ 2 + p|w||Vw| 2 K| + p\u t \ 2 \Vu\)dxdt 

9 

H\\^m\\l +||Vrf t (0)|| 2 ) + C^/, (3.17) 

i=l 
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By Lemma 13. 4[ Sobolev embedding inequality and Holder inequality, we estimate the 
term on the right hand side as follows 



h<C \\u\\i\\y^u t \\ 2 \\u t \\ 6 ds<C \\Vu\\ 2 2 \\^pu t \\ 2 \\Vu t \\ 2 ds 
Jo Jo 

<e [ WVutWlds + C [ \\Vu\\*\\^5u t \\lds 
Jo Jo 

<e f \\Vu t \\ 2 2 ds + C(\\Vu \\l + \\V 2 d \\l) 3 , 
Jo 

h<C f \\u\\l\\V 2 u\\ 2 \\u t \\ 6 ds < C [ \\Vu\\l\\V 2 u\\ 2 \\Vu t \\ 2 ds 
Jo Jo 

<e [\\Vu t \\ids + C(\\Vu \\l+ \\\V 2 d \\i)\ 
Jo 

h<C [ \\u\\l\\Vd t \\ 2 \\Vd t \\ 6 ds<C f ||V M ||2||V^|| 2 ||V 2 ^|| 2 ^ 
Jo Jo 

<e [\\V 2 d t \\ 2 2 ds + C(\\Vu \\ 2 2 + \\V 2 d \\ 2 2 ) 2 f WVdtgds, 
Jo Jo 

h<C [ \\Vd\\ 2 \\Vd\\ 6 \\u t \\lds < C [ ||W|| 2 ||V 2 d|| 2 ||V^|| 2 d S 
Jo Jo 

<CC 1/2 (||V^o|ll + II V 2 rf |ll) 1/2 f \\Vd t \\ 2 2 ds < Ce l J 2 f \\Vu t f 2 ds, 

Jo Jo 



h< \\Vd\\l\\Vd t \\ 2 \\Vd t \\,ds<C \\V 2 d\\ 2 2 \\Vd t \\ 2 \\V 2 d t \\ 2 ds 
Jo Jo 

<e f\\V 2 d t f 2 ds + C t \\V 2 d\\$\\V 2 d t \\ 2 2 ds 
Jo Jo 

<e f \\V 2 d t \\ 2 2 ds + C{\\Vu Q \\ 2 2 + \\V 2 d4l) 2 f \\Vd t \\ 2 2 ds, 
Jo Jo 

h<C [ \\u\\ 4 6 \\Vu\\ 2 6 ds < C [ \\Vu\\ 4 2 \\V 2 u\\ 2 2 ds < C(||V« ||2+ l|V 2 rf || 2 
Jo Jo 

h< f \\Vd\\t\\d t \\lds<C f \\V 2 d\\l\\Vd t \\ 2 2 ds 
Jo Jo 

<^(l|V Wo || 2 +||V 2 rf || 2 ) 2 t \\Vd t \\ 2 2 ds, 

Jo 

h<C [ \\Vpu\\ 2 \\Vu\\l\\u t \\ 6 ds < C [ \\^\\l\\V 2 u\\ 2 2 \\Vu t \\ 2 ds 
Jo Jo 

<e [ \\Vu t \\lds + C [ ||^||^||V 2 n||^ S 
Jo Jo 

<e [ t \\Vu t \\ 2 2 ds + C(\\Vu \\ 2 2 +\\V 2 d \\ 2 2 ) 3 , 
Jo 
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h <C f ||V^||2|K||e||V«||5 /2 ||V a «||J /2 tZ« 
Jo 

<C /"llV^HallVutllallVtiliy 2 \\V 2 u\& 2 ds 
Jo 

<e f \\Vutgds + C [ \\^u t \\l\\Vu\\ 2 \\V 2 u\\ 2 ds 
Jo Jo 

/ ft \V2 / ,i \ 1/2 

HVtitHlds + C / || VzzHlrfs ) / ||V 2 w|| 2 cfe sup HVp^H 2 

\Jo J \Jo J 0<s<t 



<e [\\Vu t \\lds + CCl /2 (\\Vu \\l+\\V 2 d \\l) 1/2 sup ||VM|| 2 

Jo 0<s<t 
t 



<e [ WVutgds + Ce 1 ^ 2 sup HVp^II 2 - 

Jo 0<s<t 

Substituting all these inequalities into ( 13. 17ft gives 



SUp (||^||2 : ||V«) + / '(HVTitlll + ||d*||2 + IIV^H 2 )^ 
0<s<t Jo 

<(\\VP~ou t (0)\\t + l|V*(0)||») + C(||V^o|II + \\V 2 d \\ 2 2 ) 3 + C(||Wo||l 

+ l|V 2 d ||l) 2 f \\Vd t \\lds + CC (\\Vu \\l + || V 2 d || 2 ) 2 , (3.18) 
Jo 

provided eo is small enough. Using Gagliado-Nirenberg inequality and equation (11.41) that 

||V(W| 2 d - u ■ W)|| 2 <(|||n|V 2 d|| 2 + |||V M |W|| 2 + ||Vd||jj + ||W|V 2 d||| 2 ) 

<C(||u||;|| V 2 d|| 2 || V 2 rf|| 6 + || Vd\\l\\Vu\\ 2 \\Vu\\o 

+ ||V 2 rf||^+||Vd||^||V 2 rf||2||V 2 d|| 2 ) 
<C(||V M || 2 ||V 2 rf|| 2 ||V 3 rf|| 2 + ||V 2 d||^||Vn|| 2 ||V 2 n|| 2 

+ l|V 2 rf||^||Vci|| 2 ||V 3 rf|| 2 + ||V 2 rf|| 2 ||V 2 rf|| 2 ||V 3 rf|| 2 ), (3.19) 

and thus, it follows from Lemma 13.41 that 

f \\Vd t \\lds <C A|| V 3 rf|| 2 + || V 2 w|| 2 + || Vu\\i\\V 2 d\\ 
Jo Jo 



2j||2 
2 



+ l|V 2 rf||^||V M || 2 +||V 2 rf|| 2 ||V 3 rf|| 2 )^ 
<C(l+||V M o||^ + ||V 2 rf || 2 ) 3 , 

which, combined with (13.181) . gives 

sup (II^H 2 + ||WJ 2 ) + A||VnJ 2 + ||^|| 2 + ||V 2 ^|| 2 )^ 

0<s<t Jo 

<(IIVpM0)||1 + II W £ (0)||1) + C(l + HWolll + ||V 2 d || 2 ) 4 . (3.20) 
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The compatible condition and (13.191) implies 

11^(0)1111 + ||W t (0)||l < C[\\g \\l + (1 + HVtiolll + ||Wo||^) 2 ], 
from which, recalling (I3.2ip . we arrive 

sup (WVPMI + l|V4||j|) + [\\\Vu t \\l + \\d tt \\l + \\V 2 d t \\l)ds 

0<s<t Jo 

<C[\\ga\\l + {l + \\VMl + \\VMb)% (3-21) 

By elliptic estimates to Stokes equations and elliptic equations, recalling (13.191) . it 
follows from equations (11.21) and ( j 1.4ft that 

|| V 3 d\\l <C\\Vd t \\t + e(\\V 2 u\\l + \\V 3 d\\l) + C(\\Vu\\t + || V 2 d\\t) 

and 

II V 2 u\\l + ||Vp||2 <C\\^u t \\ 2 + C(|||W|V 2 d|| 2 + \\pu\Vu\\\l) 

<C\Wph\\l + C(||Vrf||2||V 2 d|| 2 ||V 3 d|| 2 + H| 2 ||V M || 2 ||V 2 M || 2 ) 
<C||^|| 2 + e(\\V 2 u\\l + ||V 3 d|| 2 ) + C(\\V 2 d\\t + ||Vtt||S), 

and thus 

l|V 2 ^|| 2 + ||V P || 2 + ||V 3 rf|| 2 < C(\\^pu t \\ 2 2 + \\Vd t \\ 2 2 ) + C(||V 2 rf||^ + ||Vn||«). 
On account of this inequality, it follows from (13.211) and Lemma 13.41 that 
sup (UVMjj + || W t || 2 + || V 2 «|| 2 + || Vp||l + ||V 3 d|| 2 ) 

0<s<i 

+ f (l|Vu t || 2 + ||^|| 2 + ||V 2 ^|| 2 )rf S < C[\\g \\l + (1 + \\Vu \\l + ||Vd ||| 2 ) 4 ], 
Jo 

completing the proof. □ 

Proposition 3.1. Let (p,u,d,p) be a strong solution to the system U.1\) - (T77\ ) on Qt, 

and set 

E(t) = sup (||Vp|| 2 + ||pj 2 + \\u\\ 2 H2 + \\Vd\\ 2 H2 + ||Vp|| 2 ) + AllV^H 2 + ||V 2 ^|| 2 )d S . 

0<s<t JO 

Then there is a positive constant Eq depending only on ~p and Q, such that 
E(t) < C(\\ Po \\ m , |KHj*», ||Vdo||^, \\9oh,P,to) 

for any t E (0,T), provided \\^u \\l + ||Vdo||l < if N = 2, and (Wy/pHuoWl + 
\\Vd \\ 2 )(\\Vu \\ 2 2 + \\V 2 do\\ 2 2 )<e ifN = 3. 
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Proof. Using Sobolev embedding inequality and Gagliado-Nirenberg inequlity and 
applying elliptic estimates on Stokes equations, it follows from Lemma 13.41 and Lemma 
[33] that 

||V 2 «|| 2 <C{\\pu t \\l+\\pu\Vu\\\\+\\Vd\V 2 d\\\l) 

<C(\\Vu t \\l + \\Vu\\l\\ V 2 u\\l + ||V 2 d|| 2 ||V 3 d|| 2 ^ 



2) 



<C(||V^|| 2 + 1 



if JV = 2, and 



||V 2 u|| 2 <C(\\pu t \\l+ \\pu\Vu\\\l+ ||W|V 2 d||| 2 ) 

^(iiVtiti^ + lHiLIIVTiiiS+llvdi^iiv 2 ^) 

<C(||V^|| 2 + ||V 2 M || 2 ||V M ||^ 4/13 ||V 2 M ||f /13 + ||V 3 d|| 2 ||V 3 d|| 2 ) 
<^l|V 2 n|| 2 + C||Vn t || 2 + C, 

if N = 3. As a consequence, it follows from Lemma 1331 that J Q * || V 2 w||f<is < C. On account 
of this inequality, it follows from Lemma [27T1 that ||Vp(t)||| < C||Vpo||i- Using equation 
( II. ip . by Sobolev embedding inequality and Lemma 133} we deduce 

llp t ||^<^lkllLllVp|| 2 <c||v 2 M || 2 ||v P || 2 <a 

The proof is complete. □ 



4 Proof of the main theorem 

After obtaining the a priori estimates stated in the previous section, we can now give 
the proof of our main results. 

Proof. We divide the proof into several steps. Let Eq be the constant stated in 
Proposition 13.11 and suppose that 

||v^mo||^+ UVdoll! < e UN = 2, 

(llVP^oll 2 . + HVdollaXHVuoll! + || V 2 d || 2 ) < e if N = 3. 

Step 1. Global existence of strong solutions without vacuum. Suppose that there is a 
positive constant p such that po(x) > p for all x G Q. By Lemma \2A\ there is a unique 
strong solution to system (11.ip - (11.7p in Qt- Let T* be the maximal existence time. We 
claim that T* = oo. If it's not the case, then < T* < oo. By Proposition 13.11 there is a 
positive constant K, such that 

sup (llVdUi, + ||Vp|| 2 + \\u\\jp + ||Vp|| 2 + ||pi|| 2 ) + [ T (\\d t \\ 2 H2 + \\u t \\ 2 Hl )ds < K 

0<t<T Jo 

(4.1) 
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for some suitable constant K depending only on ||<?o||2, ||wo||_ff 2 ? HVdoll/p, P and Q. Since 
Po(^) > p > for all x G Q, Lemma [27T1 gives p(x,t) > p for all (x, t) G Qt,, and thus 

sup ||p~ 1/2 (-Aw-div(Vd© Vd))\\l < C 

0<t<T» 

for some constant C. Using Lemma [274| we can extend the strong solution to (0, T*) with 
some T* > T*, which contradicts to the definition of T*. This contradiction implies that 
T* = oo, and thus we obtain the global strong solutions for system (11.1j) - (11.7j) with initial 
density being away from vacuum. This completes the proof of step 1. 

Step 2. Global existence of strong solutions with vacuum. We prove it by approxima- 
tion. For any j = 1, 2, ■ • • , define p? = p + h and 

9o = (Por 1/2 (-A% - Vp - div(Vd © Wo)), 
then the compatible condition holds true 

-Au - Vj9 - div(Vrf Vd ) = \fpigl 

in Q. It's easy to see that \gl\ < \go\ in Q and thus \\gH\2 < \\g0W2 f° r & U J- Then for j 
large enough it has 

||y^Juo||2+ II V^o II 2 < ^0 if N = 2, 

(II V^oll' + l|VdoIll)(||Vw ||l + ||V 2 d || 2 ) < ^0 if N = 3. 

Using the result we have proven in step 1, for each j large enough, there is a global 
strong solution (p 7 , ,p\ dp) to system (ll.ll) - (ll.7p with initial data p^O) = p and m j (0) = 
Mo. Moreover, by Proposition 13.11 (fP , w J , p- 7 , d- 7 ) satisfies the following estimates 

sup {\\Vd\\ 2 H , + ||Vp|| 2 + Hull 2 ,, + ||Vp|| 2 + ||p t || 2 ) + f{\\d t \\ 2 H2 + < ^ 

0<t<T Jo 

By the aid of these uniform estimates with respect to j, we can follow the standard 
convergence approach to obtain the global strong solution to system (11. II) - (11. 71) satisfying 
the regularities and estimates stated in Theorem ll.il Since it's quite a standard approach, 
we omit it here. 

Step 3. Uniqueness of global strong solutions. Let (p,u,p,d) and (p,u,p,d) be two 
global strong solutions to system (11.11) to (II. 4p satisfying the regularities and estimates 
stated in Theorem ll.il We only consider the three dimensional case, the two dimensional 
case can be dealt with in the similar way. Using the identity 

(p-p)t + u- V(p - p) = {u - u) ■ Vp, 

we can deduce that 

(|p - +u ■ V(|p - pf/ 2 ) <\\u- u\\Vp\\p - pf/ 2 . 
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Integrating over Q, we deduce from Sobolev embedding inequality that 
d 



dt j Q \p- p\ v2dx AL^~ ui 1 V/}| ip ~ ^ 1/2dx 

<^\\u- uU^pUp - < C\\V{u - n)|| 2 ||Vp|| 2 ||p - 



and thus 

^IIp - p\\l/2 =^\\p - p\\l f /lj t \\p - pff/l 

<C\\V(u-u)\\ 2 \\Vp\\ 2 \\p-ph /2 

<\\\V(u-u)\\l + a(t)\\p-p\\l /2 (4.2) 

for some nonnegative function a(t) G L 1 (0,T) for any T > 0. Using equation (jl.2p . we 
deduce 

p(u — u)t + pu ■ V(m — u) — A (it — u) + V(p — p) 
=A(d - d) ■ Vd + Ad ■ V(d - d) - (p - p)(u t + u • Vu) - p(u - u) • Vu. 

Multiplying the above equation by it — w and integration by parts yields 

- — / p\u — u\ 2 dx + f |V(w — u)\ 2 dx 
2 dt Jq J a 

< [ (\V(d-d)\\Vd\\V(u-u)\ + \V(d-d)\\V 2 d\\u-u\ + |Ad||V(d - d)\\u - u\ 
Jn 

+ \p — p\ \u t + u ■ Vu\\u — u\ + p|Vw||w — u\ 2 )dx, 

consequently, it follows from Holder inequality, Sobolev embedding inequality and Young 
inequality that 

j t \\Vp(u-u)\\l+\\V(u-u)\\l 

<l\Mu-u)\\l + c\\vd\\ 2 jMd-d)g 



+ C(\\V 2 d\\ 3 \\V(d - d)\\ 2 \\u - u\\ 6 + || Ad|| 3 ||V(d - d)\\ 2 \\u - u\\ 6 
+ \\p ~ Ph/2\\u t + u ■ V«|| B ||tt - u||e + ||V-u||oo||\/p( w ~ 
<\\\V{u -u)\\l + c[{\\Vd\\lo + \\^ 2 d\\l + ||Ad||l)||V(d - d)\\l 
+ \\u t + u- Vw||^||p-p||^ /2 +||V«|| 00 ||v / p(w-«)||2 



from which we deduce 

d 



Jt \\Vp(u-u)\\i+\\V(u-u)M 
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<m(\\v(d-d)\\ 2 2 +\\ P -p\\i /2 +\\^p( 



3/2 t- ii - U)\\ 2 ), 



for some nonnegative function /3(t) G L 1 (0,T) for any T > 0. Combining the above 
inequality with (14. 2p . we obtain 



( u-u) ||i 



^(lIp-pllI/a+IIVpCti-^lllJ+HVO 

< 7 (f)(l|v(d - + Hp - p\\1/2 + - 6)111), (4-3) 

for some nonnegative function j(t) G L 1 (0, T) for any T > 0. Using (11.4)) . it has 

(d - d) t + it • V(d - d) - A(d - d) 
=V(d - d) : V(d + d)d + | Vd| 2 (d - d) - (w - u) • Vd. 



1 d 

2dty n 



Multiply the above equation by — A(d — d) and integration by parts, we obtain 

\V{d-d)\ 2 dx+ [ \A(d-d)\ 2 dx 
Jn 

= I u ■ V(d - d) + (u-u)- Vd - V(d - d) : V(d + d)d - | Vd| 2 (d - d) • A(d - d)dx 
Jn i- 

<~ / |A(d-d)| 2 dx + / [(|m| 2 + |V(d + d)| 2 )|V(d-d)| 2 + |Vd| 4 |d-d?+ |VdT| M -M| 2 
It follows from Poincare inequality and the above inequality that 

|||V(d - d)|| 2 + ||A(d - d)|| 2 < C(*)(||V(d - d)|| 2 + ||V(« - n)|| 2 ), 

for some nonnegative function ((t) G L 1 (0,T) for any T > 0. Combining the above 
inequality with (14. 3p . we have 

j t (\\p- p\\1/2 + - 6)111 + liv(d - cO Hi) + (||v(« -u)\\l+ l|A(d - 4)111) 

<e(*)(iip - p\\l/2 + - 6)in + ii v(d - ^iii), 

for some nonnegative function £(;£) G L 1 (0,T) for any T > 0. Recalling that (p — p,u — 
u,d — d)\ t =o = (0, 0, 0), it follows from Gronwall's inequality that (p, w, d) = (p, u, d), and 
consequently p = p, completing the proof. □ 
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